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We calculate the contributions arising at order a s in the QCD sum rule for the 
spin-symmetry violating universal function X3( v ' v ')i which appears at order 1/toq in 
the heavy quark expansion of meson form factors. In particular, we derive the two- 
loop perturbative contribution to the sum rule. Over the kinematic range accessible 
m B -> D^ku decays, we find that X3( v ' v ') does not exceed the level of ~ 1%, 
indicating that power corrections induced by the chromo-magnetic operator in the 
heavy quark expansion are small. 
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I. INTRODUCTION 



In the heavy quark effective theory (HQET), the hadronic matrix elements describing 
the semileptonic decays M(v) —>■ M'{v')iv, where M and M' are pseudoscalar or vector 
mesons containing a heavy quark, can be systematically expanded in inverse powers of the 
heavy quark masses |1]-§J. The coefficients in this expansion are mQ-independent, universal 
functions of the kinematic variable y = v ■ v'. These so-called Isgur-Wise form factors 
characterize the properties of the cloud of light quarks and gluons surrounding the heavy 
quarks, which act as static color sources. At leading order, a single function £ (y) suffices to 
parameterize all matrix elements ||. This is expressed in the compact trace formula |],[7| 

(MV)I J{0) \M(v)) = -ay)HM'(v')TM{v)}, (1) 

which is a consequence of a spin-flavor symmetry that QCD reveals in the limit mq — > oo 
|6|||. Here J = h'(v')Th(v) is an arbitrary heavy quark current, h(v) describes a heavy 
quark with velocity v, and 

, -p f — 75 ; pseudoscalar meson 

[y ) y/ m + | ^ • vector meson ^ ' 

is a spin wave function that describes correctly the transformation properties (under boosts 
and heavy quark spin rotations) of the meson states in the effective theory P + = |(1 + ^) is 
an on-shell projection operator, and e denotes the polarization of the vector meson. Vector 
current conservation implies that the Isgur-Wise function £ (y) is normalized to unity at zero 
recoil, where the velocities of the initial and final meson are the same, and the maximum 
momentum is transferred to the lepton pair |J. 

At order 1 / mq , the spin-flavor symmetry is explicitly broken by the presence of higher di- 
mension operators in both the effective currents and the effective Lagrangian of HQET. Some 
of these power corrections, which violate the spin symmetry, arise from the chromo-magnetic 
interaction between the spin of the heavy quark and the gluon field. The corresponding term 
in the effective Lagrangian is 

£mag = 7^ Cmag , Cmag = ^ h(v)a a pG a0 h(v) , (3) 

Iraq I 

where Z is a renormalization factor. The effects of an insertion of this operator can be 
parameterized by a tensor form factor x at3 IH3 : 



(M'(v')\ijdzT{J(0),O mag (z)}\M(v)} = -Ati{x aP (v,v')M\v')TP + ia aP M(v)}. (4) 

The mass parameter A sets the canonical scale for power corrections in HQET. In the 
mq — > oo limit, it measures the finite mass difference between a heavy meson and the 
heavy quark that it contains [|TTJ. By factoring out this parameter, x a/3 (t>,i/) becomes 
dimensionless. The most general decomposition of this form factor involves two real, scalar 
functions X2{y) and Xziy) defined by [H] 



x "P(v,v') = (</V - v'^ a ) X 2(y) - 2ia^ X3 (y). (5) 
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Irrespective of the structure of the current J, the form factor Xz{v) appears always in the 
following combination with £(y): 

^) + 2Za(^ + ^W (6) 
\rriQ rriQ'/ 

where dp = 3 for a pseudoscalar and dy = — 1 for a vector meson. It thus effectively 
renormalizes the leading Isgur-Wise function, preserving its normalization at y — 1 since 
X3(l) — according to Luke's theorem [|l(J. Eq. @ shows that knowledge of Xs(y) i s 
needed if one wants to relate processes which are connected by the spin symmetry, such as 
B -> D £ v and B -»• D*£ v. 

Being hadronic form factors, the universal functions in HQET can only be investigated 
using nonperturbative methods. QCD sum rules have become very popular for this purpose. 
They have been reformulated in the context of the effective theory and have been applied to 
the study of meson decay constants and the Isgur-Wise functions both in leading and next- 
to- leading order in the 1/toq expansion |T2|-f2T[|. In particular, it has been shown that very 



simple predictions for the spin-symmetry violating form factors are obtained when terms of 



order a s are neglected, namely [17]] 



X2(y) = o, 

xM^(99s<rapG^q)[l-^y)]. (7) 

In this approach Xziv) is proportional to the mixed quark-gluon condensate, and it was 
estimated that Xz{y) ~ 1% f° r large recoil (y ~ 1.5). In a recent work we have refined the 
prediction for Xi{y) by including contributions of order ot s in the sum rule analysis [^]. We 
found that these are as important as the contribution of the mixed condensate in (|7|). It 
is, therefore, worthwhile to include such effects also in the analysis of Xz{v)- This is the 
purpose of this article. 

II. DERIVATION OF THE SUM RULE 

The QCD sum rule analysis of the functions X2{y) and Xz{v) ls very similar. We shall, 
therefore, only briefly sketch the general procedure and refer for details to Refs. |T7lj20f . Our 
starting point is the correlator 

Jdxdx'dzeW-*'-^ (0\T{[qT M ,ti]j,J(0),O^(z),[hT M q] x }\0) 
= E 2 (u, u/, y) tr{(v' a jP - v'^ a ) T M >P+r P+ ia a/3 P + T M } 

+ E 3 {u, u', y) tx{2a aP T M ,P'+T P + a aP P + Y M ) , (8) 

where P' + = \{1 + ^')> and we omit the velocity labels in h and h! for simplicity. The heavy- 
light currents interpolate pseudoscalar or vector mesons, depending on the choice Ym = —75 
or Tm = Ifi — Vfi, respectively. The external momenta k and k' in (|||) are the "residual" 
off-shell momenta of the heavy quarks. Due to the phase redefinition of the effective heavy 
quark fields in HQET, they are related to the total momenta P and P' by k = P — rriQ v 
and k' = P' -m Q ,v' @. 
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The coefficient functions Sj are analytic in uj = 2v ■ k and uj' = 2v'-k', with discontinuities 
for positive values of these variables. They can be saturated by intermediate states which 
couple to the heavy-light currents. In particular, there is a double-pole contribution from the 
ground-state mesons M and M'. To leading order in the I/tuq expansion the pole position 
is at uj = uj' = 2 A. In the case of S2, the residue of the pole is proportional to the universal 
function Xz(y)- F° r ^3 the situation is more complicated, however, since insertions of the 
chromo-magnetic operator not only renormalize the leading Isgur-Wise function, but also 
the coupling of the heavy mesons to the interpolating heavy-light currents (i.e., the meson 
decay constants) and the physical meson masses, which define the position of the pole.Q The 



correct expression for the pole contribution to S 3 is fll? 



(u - 2A + ie)(u' -2A + ze) 

Here F is the analog of the meson decay constant in the effective theory (F ~ ju\J m M)-, G 2 
is the spin-symmetry violating correction to it, and 5A 2 denotes the spin-symmetry violating 
mass shift of the meson masses at order I/tuq. More precisely, these quantities are defined 



by [16 



m M - mn = All + — 5A 2 + . . . } 
I m Q J 



(0\j(p)\M(v)) = % -jix{rM(v)}, (10) 
(0\ij dzT{j(0), O mag (z)} \M(v)) = y G 2 tr{2^r P + a a(3 M{v)} , 

where the ellipses represent spin-symmetry conserving or higher order power corrections, 
and j = qTh(v). In terms of the vector-pseudoscalar mass splitting, the parameter 5A 2 is 
given by my — m 2 P = — 8 A 5A 2 . 

For not too small, negative values of u and uj', the coefficient function S3 can be approx- 
imated as a perturbative series in a s , supplemented by the leading power corrections in 1/u 
and l/o/, which are proportional to vacuum expectation values of local quark-gluon opera- 
tors, the so-called condensates [p2]] . This is how nonperturbative corrections are incorporated 



in this approach. The idea of QCD sum rules is to match this theoretical representation of 
S 3 to the phenomenological pole contribution given in ([]). To this end, one first writes the 
theoretical expression in terms of a double dispersion integral, 

Ef(uj, uj', y)= fdu di/ t Pf^^V) + subtractions, (11) 

J [y — uj — le) {y 1 — uj 1 — te) 

and performs a Borel transformation in uj and uj' . This yields an exponential damping 
factor, which suppresses the contamination from higher resonance states and eliminates 



There are no such additional terms for H2 because of the peculiar trace structure associated with 
this coefficient function. 
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possible subtraction terms. Because of the flavor symmetry it is natural to set the Borel 
parameters associated with uj and u' equal: r = t' = 2T. One then introduces new variables 
w± = \{v ± v') and integrates over uj_. At this stage, one employs quark-hadron duality to 
argue that the remaining integral over the "diagonal" variable uo + above a threshold Uq is 
dual to the contribution of higher resonance states fl5 , |I§ |. One thus equates the integral up 
to uq to the Borel transform of the pole contribution in (EI). This gives the QCD sum rule 



Ax 3 (z/) + 



Go 



A5A 2 
2T 



£(v) 



F 2 e 



-2A/T 



dco. 



-^/ T pl h (u; + ,y)=K(T,u; ,y). (12) 



The effective spectral density p^ h arises after integration of the double spectral density over 
uj-. Note that for each contribution to it the dependence on uj + is known on dimensional 
grounds. It thus suffices to calculate directly the Borel transform of the individual con- 
tributions to Sg h , corresponding to the limit u — > oo in fll^). The ^-dependence can be 
recovered at the end of the calculation. 

When terms of order a s are neglected, contributions to the sum rule for S 3 can only be 
proportional to condensates involving the gluon field, since there is no way to contract the 
gluon contained in O mag . The leading power correction of this type is represented by the 
diagram shown in Fig. |I|(d). It is proportional to the mixed quark-gluon condensate and, as 
shown in Ref. [17], leads to (|7]). Here we are interested in the additional contributions arising 
at order a s . They are shown in Fig. |Tj(a)-(c). Besides a two-loop perturbative contribution, 
one encounters further nonperturbative corrections proportional to the quark and the gluon 
condensate. 

Let us first present the result for the nonperturbative power corrections. We find 



K cond {T, u) ,y) 



g s (qq)T 

Qtt 



vy; J °\tJ 96vr \y + l 



48T 



(13) 



where 



r(y) 



T(n + l) 



\n(y + Jyi-l) 



dz z n e~ 



(14) 



We think it is safe to truncate the series of power corrections after the mixed condensate, 
since effects of higher dimension condensates are expected to be very small. The next-order 
power corrections would come from four-quark operators, which would contribute at the 
level of |(^)|/T 3 ~ 1 — 5% as compared to the quark condensate. 

The calculation of the perturbative contribution is more cumbersome. Two-loop dia- 
grams such as that depicted in Fig. 0(a) were considered in detail in Ref. [pT|| . We use 
dimensional regularization and the integral representations discussed in this reference. Af- 
ter Borel transformation we find 

(2T) 2D - 4 



K pert (T } oo, y) = AN C C F g 2 s T(D - V 



(4tt) 



D 



d\ X 1 - / dm / du 



l)D / 2 - 



l/A 



[ui + 2(y+l)(u 2 -l)] 



D-l ' 



(15) 
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where Cp = — 1)/2N C , and D is the dimension of space-time. For D = 4, the integrand 
diverges as A — > 0. To regulate the integral, we assume D < 2 and use a triple integration 
by parts in A to obtain an expression which can be analytically continued to the vicinity of 
D = 4. Next we set D = 4 + 2e, expand in e, write the result as an integral over u + , and 
introduce back the continuum threshold. This gives 



K pert (T, u ,y) 



o:. s 



487T 3 \y + 1 



^'0 



dcu+ u 3 e-" +/T 



e — 2 In Air + 4 In uj + — yr(y) — 2 In 



(16) 



y + 1 23 



On first sight, the appearance of a divergence at order a s seems surprising. Since the leading 
contribution to the spin-symmetry violating form factors is of order g s , one would expect di- 
vergences to appear at order g 3 . In fact, the one-loop renormalization group equations which 
control the running of the hadronic quantities Xs(y) an d SA 2 are homogeneous equations, 



confirming this assertion. But G 2 satisfies the inhomogeneous equation (JT6 

d 



^A. 
9tt 



(17) 



which shows that divergences arise at order a s . At this order, the renormalization of the 
sum rule is thus accomplished by a renormalization of the "bare" parameter G 2 in (|12"D. In 
the MS subtraction scheme, one defines a renormalized parameter G 2 (fi) by 



G 2 {fi) 



,^ybare , . 

° 2 + 9^ A 



1 



+ 7b - In 



47T\ 



0{gl) . 



(18) 



Hence a counterterm proportional to A£(y) has to be added to the bracket on the left-hand 
side of the sum rule (12). To evaluate its effect on the right-hand side, we note that in D 



dimensions 1 7 



M{y) F 2 e~ 2R / T 



lQ7i 2 \y + l 



(19) 



x 1 + e 



7£ — In 47T + 2 In uj + — In 



y + 1 



- 2 



} + 0(g s ,e 2 ) 



From fll2]) , (|16D, (p~8[ ) and (|T9"|) it is seen that indeed the 1/e pole cancels upon renormalization 
of G 2 . We find 



Ax 3 (z/) + 



G 2 (/x) 



A5A 2 



2T 



tt.s 



48tt 3 + 1 



+-^cond(T, w ,Z/) 



-u>+/T 



, A 4 17 
2 In — + — 
6 



- 1 + in 



y + i 



(20) 



In this expression, all hadronic parameters are defined at the scale /1. However, we have 
only made the /z-dependence of G 2 {^j) explicit since this is what we are sensitive to in our 
analysis. 
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According to Luke's theorem, the universal function Xsiu) vanishes at zero recoil [10 . 



Evaluating ( P0j) for y = 1, we thus obtain a sum rule for G 2 (/i) and <5A 2 . It reads 



C 2 (m) - 



A5A 2 
2T 



^2 e -2A/T 



a. 



24vT 3 



hi 
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+ K coad (T,uj ,l), (21) 



where we have used that r(l) = 1. Precisely this sum rule has been derived previously, 
starting from a two-current correlator, in Ref. |Tj|. This provides a nontrivial check of our 
calculation. Using the fact that £(y) = [2/(y + l)] 2 + 0(g s ) according to (0), we find that 
the /^-dependent terms cancel out when we eliminate G 2 (/U and 5A 2 from the sum rule for 

Before we present our final result, there is one more effect which has to be taken into 
account, namely a spin-symmetry violating correction to the continuum threshold ujq. Since 
the chromo- magnetic interaction changes the masses of the ground-state mesons [cf. ([10|)1, 
it also changes the masses of higher resonance states. Expanding the physical threshold as 



^'pi.vs = uj \ 1 + —5uj 2 + 



(22) 



we expect 8u 2 — 5A 2 — —0.12 GeV, where the numerical value follows from the observed 
mass splitting between B* and B. The appearance of such a shift is also in the spirit 
of QCD sum rules. Assume that we had constructed a sum rule for the combination of 
form factors given in (|6]). Then the presence of 1/tuq corrections would affect the stability 
of the sum rule for the Isgur-Wise function £(y) alone. The parameters that guaranteed 
optimal stability for £(?/) would have to be readjusted by an amount of order 1/mq in order 
to provide optimal stability of the sum rule for the combination @. From such a self- 
consistent analysis, the parameters <5A 2 and <5u; 2 have been determined in Ref. [IB|. One 
finds the vector-pseudoscalar mass splitting in excellent agreement with experiment, and 
5u>2 = —(0.10 ± 0.02) GeV, which is just what we expected on physical grounds. 

The contribution of 8u 2 to the sum rule for Xziy) has been calculated in Ref. 
Including this term, we now present our final result: 
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Xs(y)AF 2 e- 2A / T 



a s T 4 



y + 1 



.^2 3 uo/T 
.2 ^0 e 



yr(y) — 1 + In 
2 ^ 2 



y + 1 



32tt 

6n 
(a s GG) 



96vr 



y + 1 

± -r{y)-i{y) 
2 



2 



y + 1 



ay) 



(qgs<7af>G a Pq) 
48T 



(23) 



It explicitly exhibits the fact that X3(l) — 0. 



III. NUMERICAL ANALYSIS 

Let us now turn to the evaluation of the sum rule (|23|) . For the QCD parameters we 
take the standard values 
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(qq) = -(0.23 GeV) 3 , 
(a 8 GG) = 0.04 GeV 4 , 
(qg s a a pG a0 q) = m 2 (qq) , m 2 = 0.8 GeV 2 . (24) 

Furthermore, we use 5u)2 = —0.1 GeV from above, and a s /ir = 0.1 corresponding to the scale 
fx = 2A ~ 1 GeV, which is appropriate for evaluating radiative corrections in the effective 
theory fll5|[ . The sensitivity of our results to changes in these parameters will be discussed 
below. The dependence of the left-hand side of ( ^3|) on A and F can be eliminated by using 



a QCD sum rule for these parameters, too. It reads [16 



t F 2 -2A/T = 9T^ r ("0 \ _ (qg s CT a pG a Pq) 
" 8tt 2 3 VT/ AT 



(25) 



Similarly, we calculate the Isgur-Wise function from £(y) = F(T, u , y)/F(T, u , 1), where 

El 

-ptrp s 3T 3 / 2 y f u , (2y + l) (gg s a af3 G^g) 

-o, y) = — 2 5 2 (-) - (qq) + — — . (26) 



Combining (|23|) , (g5|) and (pq) , we obtain X3(y) as a function of a;o and T. These parameters 
can be determined from the analysis of a QCD sum rule for the correlator of two heavy-light 
currents in the effective theory p!6|JT8| . One finds good stability for ujq = 2.0 ±0.3 GeV, and 
the consistency of the theoretical calculation requires that the Borel parameter be in the 
range 0.6 < T < 1.0 GeV. It supports the self-consistency of the approach that, as shown 
in Fig. H], we find stability of the sum rule (|23|) in the same region of parameter space. Note 
that it is in fact the (^-term that stabilizes the sum rule. Without it there were no plateau. 

Over the kinematic range accessible in semileptonic B — > D^*'iv decays, we show in 
Fig. ^|(a) the range of predictions for Xz{v) obtained for 1.7 < u < 2.3 GeV and 0.7 < T < 
1.2 GeV. From this we estimate a relative uncertainty of ~ ±25%, which is mainly due to 
the uncertainty in the continuum threshold. It is apparent that the form factor is small, not 
exceeding the level of 

Finally, we show in Fig. |3](b) the contributions of the individual terms in the sum rule 
(f2~3|). Due to the large negative contribution proportional to the quark condensate, the terms 
of order a s , which we have calculated in this paper, cancel each other to a large extent. As 
a consequence, our final result for Xziv) ls n °t very different from that obtained neglecting 



these terms ||17|| . This is, however, an accident. For instance, the order-a s corrections would 
enhance the sum rule prediction by a factor of two if the (gg)-term had the opposite sign. 
From this figure one can also deduce how changes in the values of the vacuum condensates 
would affect the numerical results. As long as one stays within the standard limits, the 
sensitivity to such changes is in fact rather small. For instance, working with the larger 
value (qq) = —(0.26 GeV) 3 , or varying m 2 , between 0.6 and 1.0 GeV 2 , changes Xsiu) by no 
more than ±0.15%. 



2 When comparing our result to the function X3 6n (y) shown in Fig. 2 of Ref. [O], one has to include 

— 1/3 

a renormalization factor, which is approximately given by a s — 1.5. 
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In conclusion, we have presented the complete order-a s QCD sum rule analysis of the 
subleading Isgur-Wise function Xz{y)i including in particular the two-loop perturbative con- 
tribution. We find that over the kinematic region accessible in semileptonic B decays this 
form factor is small, typically of the order of 1%. When combined with our previous analysis 
20| , which predicted similarly small values for the universal function Xi{y)i these results 



strongly indicate that power corrections in the heavy quark expansion which are induced 
by the chromo-magnetic interaction between the gluon field and the heavy quark spin are 
small. 
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FIGURES 



FIG. 1. Diagrams contributing to the sum rule for the universal form factor X2,{ v ' v ') : two- loop 
perturbative contribution (a), and nonperturbative contributions proportional to the quark con- 
densate (b), the gluon condensate (c), and the mixed condensate (d). Heavy quark propagators 
are drawn as double lines. The square represents the chromo-magnetic operator. 

FIG. 2. Analysis of the stability region for the sum rule (^): The form factor Xs(y) ls 
shown for y = 1.5 as a function of the Borel parameter. From top to bottom, the solid curves 
refer to ujq = 1.7, 2.0, and 2.3 GeV. The dashes lines are obtained by neglecting the contribution 
proportional to 8u)2 ■ 

FIG. 3. (a) Prediction for the form factor X3( v " v> ) m the stability region 1.7 < ujq < 2.3 
GeV and 0.7 < T < 1.2 GeV. (b) Individual contributions to xziy ■ v') for T = 0.8 GeV and 
u>o = 2.0 GeV: total (solid), mixed condensate (dashed-dotted), gluon condensate (wide dots), 
quark condensate (dashes). The perturbative contribution and the <5cj2-term are indistinguishable 
in this figure and are both represented by the narrow dots. 
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